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Abstract 

Even after the negative result by the MiniBooNE experiment, we can still have a (3+l)-scheme 
with one sterile neutrino whose mixing lies within the allowed region by MiniBooNE. In this note 
we discuss the possible effects of such a scheme on the flavor ratio of high-energy cosmic neutrinos 
from cosmologically distant astrophysical sources. It is shown that in principle there is still a 
chance to observe deviation from the standard three flavor scenario in the flavor ratio of the high- 
energy cosmic neutrinos. It is proposed to see the energy spectrum to cope with the theoretical 
uncertainties which were recently pointed out by Lipari, Lusignoli and Meloni. It is emphasized 
that the v„ — u T ratio is relatively insensitive to the theoretical uncertainties and therefore this 
ratio is the key ingredient to look for the signatures of the sterile neutrino scheme. Although 
the statistics of data from one source in the next generation of neutrino telescopes are estimated 
not to be sufficient to distinguish the three and four family schemes, if we can gain statistics by, 
e.g., summing over data from many sources, then it might be possible to have a signature for the 
(3+l)-scheme. 
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I. INTRODUCTION 



There has recently been much interest in high energy cosmic neutrinos from cosmolog- 
ically distant astrophysical sources, e.g., Active Galactic Nuclei (AGN) and Gamma Ray 
Bursts (GRB) (see, e.g., motivated by their possible observation at ongoing or future 
experiments such aslGECUBE Q, BAIKAL [3|, NESTOR [4], ANTARES [4 NEMO fl, 
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flux has been extensively studied by many people 

since the possibility of flavor tagging of high energy neutrinos was advanced in Ref. 
A measurement of the flavor composition can, in principle, give important informations on 
the properties of the astrophysical neutrino source, on the leptonic mixing sector or point 
out new physics phenomena beyond the Standard Model of fundamental interactions. This 
claim, however, has been critically examined in Ref. [22), where it was shown that it will 
be difficult to learn something on the neutrino sources with the expected precision on high 
energy cosmic neutrino flavor measurements at ongoing and planned experiments. To illus- 
trate the point, the authors of Ref. (22I consider in detail the case of neutrinos produced by 
7 ray bursts (GRB) and take into account their energy distribution according to Waxman 
and Bahcall 23j ( in the framework of the "fireball model" of GRB, see for example Ref. 24] 
and Refs. therein). The outcome of their analysis is that it is not clear if discussions on the 
astrophysical neutrino flavor composition make sense at all. 

We do not pretend in this paper to clarify if a measurement of the neutrino flavor fluxes 
is possible in the absence of a precise knowledge of the source. We, more modestly, present 
an update of the analysis that was performed in Ref. 12J on the possibility to use such 
measurements at neutrino telescopes to constrain four neutrino models. 

Historically, four neutrino mass schemes became popular after the LSND group an- 
nounced their data which suggest neutrino oscillation with the mass squared difference 
Am 2 



0(1) eV 2 
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271 ] . The fourth neutrino has to be sterile, because the number of 

. In Ref. Q, the 



the light active neutrinos is three to be consistent with the LEP data 
neutrino telescope flavor data were used to constrain one specific four neutrino model, the 
so-called (2+2)-scheme, that could explain the LSND data together with negative results 
at other experiments. Recently, the MiniBooNE experiment [29j gave a negative result for 
neutrino oscillations with the mass squared difference Am 2 ~ 0(1 ) eV 2 which was suggested 



by the LSND data. While several scenarios have been proposed [30j, [31|, [32j, [33j to account 
for both the affirmative result by LSND and the negative one by MiniBooNE, four neutrino 
models are currently unable to explain the LSND data. 

However, sterile neutrino scenarios which satisfy all the experimental constraints except 
LSND are clearly possible. In Ref. [34J, a detailed analysis was performed on the allowed 
region of the mixing angles of the so-called (3+l)-scheme, which satisfies all the constraints 
but LSND. In this paper we perform an analysis of the observed ratio of the high energy 
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neutrinos in the case of the (3+l)-scheme for the allowed region shown in Ref. [34(, and 
discuss whether we can establish the signature of this scenario by looking at the ratios 
Reix = {ye + Ve)/{vn + ^v) an d R T [i = {y T + v T )/{v^ + Vy) of the neutrino flux taking into 



account the uncertainties discussed in Ref. [22 ] . 



A similar analysis has been carried out in Ref. 19) . The main differences between this 
paper and Ref. [l9|] are the following: 

• We leave all the parameters of the (3+1) four neutrino model (both angles and CP- 
violating phases) free to vary, but only within the presently allowed region, as from 



Ref. [34|. We point out the critical role played by #34 in the difference between three- 
and four-family mixing. However, since #34 cannot exceed ~ 35°, no striking sterile 
neutrino signal has to be foreseen at a neutrino telescope experiment. Searching sterile 
neutrinos at this generation of experiments will be a painful task, if possible at all. 

• We study the energy dependence of the neutrino flavor fluxes and of the flavor ratios in 
the case of neutrinos from GRB, following the analysis of Ref. [zj. We show that the 
uncertainties induced by the neutrino flux energy spectrum in R efl makes impossible to 
use this observable to distinguish three- from four-family mixing. On the other hand, 
we find that R TfJ , is much less affected by the energy dependence. The flavor ratio R Tfi 
is the optimal observable to look for sterile neutrinos. 

• We present analytic formulae for the neutrino flavor fluxes and ratios expanded at 
second order in the small parameters of the model, #13, #14, #24 and 77 = 7r/4 — 2 3- 
These formulae help in the understanding of the numerical results presented in the 
paper. 

The paper is organized as follows. In Sect. [IT] we briefly summarize the (3+l)-scheme 
and the allowed region of the mixing angles which are given by all the experimental data 
except LSND. In Sect. [ITT] we study the effects of the (3+l)-scheme on the flavor ratios 
Ren — (y e + v^jiy^ + Vy) and R TfJi = {u T + u T )/{u^ + v^). In Sect. |TV]we discuss the energy 
dependence of the flavor ratios in the case of neutrinos produced in GRB's and examine 
whether deviations from the three-family expected values are larg e enough compared with 
the theoretical uncertainties of the neutrino flux given in Ref. 22|. In Sect. |V]we present a 
short discussion on the statistical requirements needed to identify a sterile neutrino signal 
using the flavor ratios. In Sect. I VI] we eventually draw our conclusions. Analytic formulae 
at second order in # 13 , #14, #24, "T? = tt/4 — #23 and in the energy-dependent parameter e(E v ) 
that is introduced to describe possible deviations from (u e , z/ T ) = (1,2 — e, 0) are shown 
in App. [A] 
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II. FOUR NEUTRINO MASS SCHEMES 



Four-neutrino schemes are classified into two classes: (3+1)- and (2+2)- schemes, de- 
pending whether one or two mass state (s) are separated from the other by a 0(1) eV 2 
mass-squared gap. When one tries to account for the LSND data and for the negative re- 
sults at atmospheric and short baseline experiments with four-neutrino schemes, it has been 
found that we fall into problems irrespectively of which scheme we choose. 

In (2+2)-schemes, the sterile neutrino contributes either to the solar or atmospheric os- 
cillations. The fraction of sterile neutrino contributions to solar and atmospheric oscillations 
is given by \U s i\ 2 + \U S 2\ 2 and \U s3 \ 2 + |t/ S 4| 2 , respectively, where the mass squared differ- 
ences Amjj and |Am! 3 | are assumed to be those of the solar and atmospheric oscillations. 
The experimental results show that mixing among active neutrinos give dominant contri- 
butions to both the solar and atmospheric oscillations (see, e.g., Ref. 35|). In particular, 
in Fig. 19 of Ref. Q we can see that at the 99% level r] s = \U sl \ 2 + \U s2 \ 2 < 0.25 and 



u s ,\ 2 + \u : 



.s i I 



sj I 



1. 



2 < 0.25, which contradicts the unitarity condition Y2j=i 
In fact the (2+2)-schemes are excluded at 5.1a CL 35j. This conclusion is independent of 
whether we take the LSND data into consideration or not. 

In the (3+l)-scheme, on the other hand, in order to account for LSND without contradic- 
tion with other disappearance experiments, the oscillation probabilities of the appearance 
and disappearance channels have to satisfy the following relation (361 . 37]: 



sin 2 2^ LSND (Am 2 



< - sin 26 l Bugcy(Am ) ■ sin 2# CDHSW (Am 



(1) 



where 6 l LSND (Am 2 ), 6 l CDHSW (Am 2 ), 6 , B Ug ey(Am 2 ) are the value of the two flavor mixing angle 
as a function of the mass squared difference Am 2 in the allowed region for LSND [p l 



and the Bugey experiment [3 




From 
gives 



the CDHSW experiment [38 1 

the early stages of the four neutrino schemes it has been known 
a strong constraint because there is very little region of the value of Am 2 which satisfies 
Eq. ([!]). While the significance to exclude the (3+l)-scheme changed 40j when the LSND 
allowed region shifted slightly toward a smaller mixing region, difficulty to satisfy Eq. (fTJ) 
is basically the reason why this scheme is also disfavored. 

A (3+2)-scheme with two sterile neutrino has also been proposed |4l| to account for 
LSND, and it may be possible to reconcile the LSND and MiniBooNE data by introducing 



a CP phase 42|, |43] . However, tension with the disappearance experiments such as CDHSW 



and Bugey always remains, as long as we take into account the LSND data. 

In this paper we will consider a (3+l)-scheme without taking the LSND data into account, 
e.g., we will assume a (3+l)-scheme which satisfies all the negative constraints given by the 
appearance experiments as it was done in Ref. 34|. Then we no longer have the constraint 
(TT]) and we have only the upper bound on the extra mixing angles, as we will see below. 
Apart from our strategy not to take the LSND constraint into account, our framework is 
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the same as the conventional (3+l)-scheme. 

There has been discussions in which the mixing angles of four neutrino schemes may be 



constrained by big-bang nucleosynthesis [3g,|44], and if such arguments are applied, then the 



mixing angles of sterile neutrinos would be very small, and no significant signature would 
be expected in the flavor ratio of high energy cosmic neutrinos. However, it was shown in 
some model 45j that neutrino oscillations themselves create large lepton asymmetries which 



prevent sterile neutrinos from being in thermal equilibrium, so it is not so clear whether 



the arguments in 3g, |44j hold. At present, therefore, it is fare to say that there is not yet 



general consensus on this issue (See Ref. [46j| and references therein.). In this paper we will 
not impose cosmological constraints on our scheme. 

In four neutrino schemes, the 4x4 mixing matrix U contains six mixing angles 6{j 
and three Dirac CP phases 5j (where we ignore the three Majorana CP phases, since only 
Dirac CP phases can be measured in oscillation experiments). There are plenty of different 
parametrizations of a general unitary 4x4 mixing matrix U. The following parametrization 
[43( | . valid in the "atmospheric regime" (i.e. where Am% ol L/E — > 0, Am1 tm L/E ~ 0(1) and 
\Am^ BL L/E\ 3> 1), is particularly helpful to extract the allowed region of the parameters of 
the present (3+l)-scheme from experiments: 

U = #34^34) #24^24) #23(023 ; £3) #14(014) #13 (#13 j #12(012, #l) • (2) 

In Eq. ([2D R jk (6) = exp[9A jk ] and R jk (9,8) = esxp[-i{S/2)D jk ]exp[9A jk ]exp[i(5/2)D jk ] 
are 4 x 4 rotation matrices, where i(A jk ) im = i{8 jt - 5 km ) and {D jk ) tm = 5 je 5 im - S M 5 im 
are the generators of the stt(4) Lie algebra. With this parametrization the phase 8% drops 
in the limit Am| x — > 0, U reduces to the mixing matrix in the standard three flavor case 
when 614, 824, 034 — > and all the phases disappear from the oscillation probabilities in the 
one-mass dominance limit. 

The mixing matrix elements in this parametrization are: 





C12C13C14 

Cl3Cl4Sl 2 e" 

cus l3 e- t52 



- S14 

-C23C 2 4Sl 2 e i<51 - C12 [c24S 13 S2 3 e l{S2 ~ S3) + C13S14S24] 
C12C23C24 - S 12 e~ lSl [c 2 4Si3S23e j(52_<53) + C13S14S24] 



ci 3 c 2 4S23e 163 - S13.S14.S24e 162 
C14S24 



(3) 



(4) 
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U t3 



Si 2 e iSl [c 34 s 2 3e i<53 + P23S24S34] 

-Cl2 {C13C24S14S34 + S 13 e l52 [C23C34 - S23S24S 3 4e" J<53 ] } 
-Cl2 [c 34 S 2 3e J<53 + C23S24S34] 

-si 2 e" i ' 51 {C13C24S14S34 + si 3 e i(52 [C23C34 - s 23 S24S34e"* <53 ] } 



(5) 



-i$2 



+ C13 [C23C34 - S 23 S24S 3 4e th ] 



U, 



«2 



where c,- 



cos 6j 



-C 2 4Sl 3 Sl4S34e 
C14C24S34 

si 2 e i<51 [C23C34S24 - s 2 3S34e j53 l 

-C12 {C13C24C34S14 - Si 3 e 4<52 [c 3 4S23S24e"*' 53 + C23S34] } 
-Cl2 [C23C34S24 - S23S 3 4e J<53 ] 

-s 12 e~ lSl {C13C24C34S14 - _s 13 e l52 [c 3 4S23S24e" 1 ' 53 + C23S34] } 

— C24C34Si3Si4e 

C14C24C34 



(6) 



C13 [c 34 S23S24e 163 + C23S34] 



and Si 



sin 6j 



The constraints from the affirmative data of the solar, KamLAND and atmospheric neu- 
trino experiments, and the negative results of the short baseline experiments on the (3+1)- 
scheme were examined in detail in Ref. [34( and the allowed region for #i 3 , 9u, 624, #34 are 
given in Fig.2 of Ref. js^]. From this figure we can see that 6*13, 6*14 and 6*24 are constrained 
to be smaller than ~ 10°, whereas #34 can be as large as 35°. A fourth "small" parameter 
that will be useful for the analytic understanding of our results is rj — 7r/4 — #23, constrained 
by atmospheric and LBL data to be \rj\ < 10°. In the analytic formulae we will consider 
#13 j #14, #24 and 77 being of the same order and expand in powers of the four. 



III. THE OBSERVED FLUXES OF EACH FLAVOR 



The oscillation lengths of standard neutrino flavor transitions are rather short with respect 
to the typical astrophysical distances. It is, therefore, a good approximation to average over 
the oscillation frequencies. The oscillation probabilities, thus, have an extremely simple 
expression, 

TV 

P a p = P{u a ^u p )=^2\U ai \ 2 \Up i \ 2 , (7) 

i=i 

where N stands for the number of neutrino flavors. It can be seen that P a p does not depend 
on the neutrino energy nor on the distance from the source. The neutrino flux of flavor a 
at the detector is: 

/3=e,/i,T 

where <3?° is the neutrino flux of flavor a at the source. Notice that, in the absence of a 
mechanism to produce sterile neutrinos at the source, the sum runs only over active neutri- 
nos. However, in the four neutrino scheme discussed later, we expect a sterile component at 
the detector (albeit small due to the severe bounds on the active-sterile mixing angles). 
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In the standard picture of neutrino creation from pion decays after interactions between 
accelerated protons and photons, the flavor composition at the source is expected to be 



approximately $° : $° : $° = 1 : 2 : 0, as in the case of the atmospheric neutrinos. If this 



' A' 



argument holds, then the neutrino flavor fluxes at the detector are, therefore, 



N 



= l^l 2 W + 2 \ U ^] $ e , (9) 



i=l 



with $° the flux of electron neutrinos at the source. 

In the case of the standard three flavor mixing (N = 3), when #13 is small and the 
atmospheric neutrino mixing angle is almost maximal, |0 2 3 — 7r/4| -C 1, it is known [10] that 
the ratio of the flavor fluxes is approximately given by 

$ e : $ M : $ r ~ 1 : 1 : 1. (10) 



In the triangle representation of neutrino flavor fluxes, proposed in Ref . [12J , the distance of 
a point from each of the sides of a regular triangle represents the fraction of neutrinos of a 
given flavor, $ Q , to the total neutrino flux, Yl a The allowed region in this representation 
for the three-family case is given by the red thin solid line in Fig. [IJ where the allowed region 



at 90%CL in Ref. [35j was used. 

In the case of four neutrino schemes (N = 4) the plot should in principle be three 
dimensional, because the flux of sterile neutrinos has to be also taken into consideration. 
As far as experiments are concerned, however, it is only the three active neutrinos that we 
can observe. We normalize, thus, the flux of each active neutrino by the total one of active 
neutrinos, so that we can still use the triangle representation also in the four flavor case. 

In 2000, when the (2+2)-scheme as well as all the solar neutrino solutions were still 
acceptable, the allowed regions in the triangle representation were relatively large (See Fig. 4 
in Ref. [l2|). One of the goal of this paper is to update that analysis to the case of the 
(3+l)-scheme, when all experimental data but the LSND ones are taken into account. We 
have performed both numerical and analytical studies to obtain the allowed region in the 
triangle representation for the (3+l)-scheme, when all the mixing angles #13, #14, #24 and 



834 are left free to vary within the allowed region in Fig. 2 of Ref. 34]. 

The allowed regions at 90% CL (blue thick solid line) and at 3a CL (green dashed line) 
for the (3+l)-scheme are shown in Fig. [TJ where the numerical results have been obtained 
assuming the standard ratio of initial neutrino fluxes, $° : $° : $° = (1 : 2 : 0). As it 
can be seen, deviations of the (3+l)-scheme flavor fluxes from the three-family case are still 
noticeable 1 . 



In Ref. [221 ] . it has been pointed out that for a neutrino flux produced by pion decay, the correct flavor 
ratios are 3?° : $° : <I>° = (1 : f .86 : 0). We have numerically checked that the allowed regions in this case 
are practically identical to those shown for the (1:2:0) case. 
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e 




FIG. 1: The allowed region of the flavor ratio <l> e : ^ : <5 T in the triangle representation 12] for 
the three flavor case (bounded by the red thin solid line) and the (3+l)-scheme (bounded by the 
blue thick solid line) at 90% CL. The standard ratio of the initial flux $° : <3?° : $° = 1 : 2 : is 
assumed. For the four flavor case the region at 3a CL (bounded by the green dashed line) is also 
depicted. The black dot in the center indicates the case for $ e : : $ T = 1 : 1 : 1. 

The results shown in Fig. [1] are easily understood using the results from our analytical 
study. At first order in the small parameters 6*13, #14, 624 and r\ = 7r/4 — 6*23, we get: 

$° [1 + (77 + S13 cos if cot 20i 2 ) sin 2 26 12 ] 
1 - § (77 + S13 cos (p cot 20 12 ) sin 2 26 12 ] 




e 
e 

r 2 ■ 



where = <5i — 62 + £3 is the only combination of CP-violating phases present in this 
approximation. Clearly, the sum of the four fluxes gives ^ Q $ Q = 3$°, as it should be in 
the absence of neutrino decays. We can see that, at first order, deviations from (1:1:1:0) are 
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governed by #13,77 and #34. The only effect of active-sterile mixing arises, thus, through #34, 
on which the experimental bounds are the least severe (#34 can be as large as 35°). A large 
#34 can make a non trivial difference between $ M and <J> T that allows us to distinguish the 
(3+ f) scheme from the three-family case. If, on the other hand, #34 were as small as #14 and 
#24, at leading order we would get a "three-family" -like result, (z/ e , v Ty v s ) = (1 + 0(e) : 
1 + 0(e) : 1 + 0(e) : 0(e 2 )). Eventually, in the limit # 34 -> 0, Eq. §TQ is reduced to the 
three flavor case with (p = 5. 

Expressions for the flavor fluxes at second order in the small parameters #13 , #14 , #24 and 
77 are shown in Appendix |A] Here, we only point out that a dependence of the fluxes on 
#14, #24 and CP- violating phases other than <p arises at this order. 



IV. ENERGY DEPENDENCE OF THE FLAVOR FLUXES 



Recently, the authors of Ref. 22| re-examined critically the theoretical uncertainties in 
the predictions for the flavor ratio as well as its energy dependence. It has been shown that, 
quite generally, an energy dependence is induced in the neutrino flavor fluxes by the specific 
characteristics of the source. To illustrate their point, a detailed discussion of the energy 
dependence of neutrino flavor fluxes originated from GRB is given. We will follow in this 
section the approach of Ref. j^] to determine the energy dependence of the flavor ratios 
in presence of sterile neutrinos. 

The analysis is performed following Waxman and Bahcall (see Ref. [23| ) for neutrinos 
from GRB. The GRB itself is described in the so-called "fireball model" [24|. Notice that 
an alternative theory has been proposed to explain 7-ray bursts, the so-called "cannonball 



model" [47| . In Ref. [481 ]. however, it was shown that no significant neutrino flux is to be 
expected from GRB's in the framework of the cannonball model. 

Following Ref. (23]], the neutrino flux is supposed to mimic the photon flux 2 and to 
have a broken power-law spectrum. The neutrino flux depends on several phenomenological 
parameters that are needed to describe the (unknown) properties of the source: (i) The 
proton spectrum N p (E p ) is assumed to depend on the proton energy E p : 



N P {E P ) oc E 



a 
V ' 



(ii) The energy dependence of the number density n 7 (e) of the photons in the wind is 
described by the two exponents of the photon energy e: 

(e/e b )"ft for e < e b , 



n 7 (e) oc I (e/e b )-^ 2 for e b < e < 

c max, 



for e > e 



max ■ 



2 We will not enter in the debate concerning the validity of this hypothesis, though, since we only use the 
neutrino spectrum obtained by Waxman and Bahcall to show how the energy dependence of the neutrino 
flavor ratios makes extremely difficult to distinguish new physics effects from standard three-family mixing. 
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where eb is the break energy and e max is the cutoff energy, (iii) The energy loss of muons 
due to synchrotron radiation is characterized by the magnetic field B which is expressed in 
terms of the parameter e M 



^ = 8.4 x 10 4 ( ^) (J*-) 
E* \ B J \KeVJ 



where E* ~ 6.9 x 10 13 (eb/KeV) _1 eV is the proton threshold energy for inelastic interactions 
with photons having the break energy eb, and e M = oo corresponds to the case with no 
energy loss. For e b ~ 1 KeV, we have E* = O(10 5 ) GeV. 

When the parameters (a, Pi, /3 2 ) are varied around the standard values (a — 2, (3\ — 1, 
/?2 = 2), the ratio ($g + $g)/($° + at the source becomes energy dependent, where <3>^ 
stands for the flux of antineutrino of flavor a at the source, and the ratio has significant 
deviation from 1/2 for lower and higher values of the neutrino energy (cf. Fig. 14 in Ref. 



22J). When the muon energy loss due to synchrotron radiation is turned on, an additional 
dam ping effect occurs and the ratio R® decreases at high neutrino energy (cf. Fig. 15 in 



Ref. |22|). 

We have studied the impact of these effects on the flux ratio which is observed on the 
Earth in the case of the three and four flavor mixing schemes. In general, if the flux ratio 
at the source is: 

(*2 + *!!) = K + $ ?) : ( $ ° = ! = KEu) : 0, (12) 

where A is a function of the neutrino energy defined as 

$° + $9 



M 1 V 



$0 + $0 > 

then the flavor ratios observed on the Earth can be written as 

$ e + $g Pee + AP^ e 

eM ~ $ +$- ~ P +XP 1 j 

^ fi ~ ^ ft 1 efi ~ /x± lift 
(J) _|_ p _i_ \ p 

TM ~ $ +$- ~ P +AP { 1 

^ ft ~ ^ ft 1 efi ~ nl fifi 

where we have used the properties P a/ 3 = Pp a and P Sj g = P{p a — > vp) = J2j Waj\ 2 Wi3j\ 2 = 
P a /3 which can be easily derived from Eq. (jTJ). 

Although the flux ratios of v e j v e and v^j are in general different from unity, the only 
necessary information to compute the flux ratio observed on the Earth is the value of A. The 
value of 1/A = R\ is plotted in Figs. 14 and 15 in Ref. [22j as a function of the neutrino 
energy for the two cases mentioned above (i.e. without and with the muon energy loss due 
to synchrotron radiation, respectively). 

From those figures, we can see that the energy- dependence of the muon-to-electron ratio 
at the source is rather small for most of the considered energy range. We can thus introduce 
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FIG. 2: The observed ratios {y e + v^jiy^ + v^) and + P T )/(z/ M + z^) for neutrinos from a GRB 
source as functions of the neutrino energy. The five cases (the three flavor case with the best- 
fit oscillation parameters; the three flavor cases which give the maximum and minimum ratios; 
the four flavor cases which give the maximum and minimum ratios) are considered with the four 
reference parameters (a = 2,j3\ = 1,02 = 2) ; (a = 2,(3% = 0.65, 02 = 2), (a = 2, fli = 1, (3 2 = 2.4), 
(a = 2.4, /?i = l,/32 = 2) in Fig. 14 of Ref. 



22 1 for each case. 



a new parameter, e(E v ) = 2 — X(E U ), and expand in #13, #14, #24, 77 and e. At first order in 
these parameters, we get from eqs. (fl~3l) and ([1] 



p(4-fam) 
1 L efj, 



1 + 2 U 



^ sin 2 29 12 ^j + ^ (<7 + On cos <p cot 12 ) sin 2 2# 12 = i?gr fam) (15) 

(16) 



p(4-fam) _ 2 _ 2 p(3-fam) 
-""r/x ~~ c 34 ~~ t ^ JX TH 

Expressions for R e ^ and _R rAt at second order in # 13 , # 14 , 6*24, 77 and e can be found in App. DAI 
The relation between three- and four-family ratios is immediately understood if we notice 
that, at first order in small parameters, we have the following simple relations between the 
mixing matrix elements |f/ a? | 2 for four flavors and |C/5 | 2 for three flavors: 
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(J 


= 4) 


r 2 \U {3) 
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(J = 1,2,3) 


s 2 

b 34 




(j = 4) 


s 2 ir/ 3) 

b 34:\ U Tj 


2 


(j = 1,2,3) 


r 2 
c 34 




(J = 4) 



(17) 
(18) 
(19) 
(20) 



where the combination of the CP phases ip = 5\ — 62 + £3 becomes 5 in the standard 
parametrization in the three flavor case. 
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From eqs. (fT5l) and (TT6j) we see that an energy dependence in the Ve/v^ ratio arises 
already at first order in small parameters, whereas the Vt/v^ ratio is less sensitive to the 
energy dependence of the source. On the other hand, deviations from the three-family 
scenario should be bigger in the latter ratio than in the former, as a consequence of the less 
stringent bounds that we have on # 34 with respect to #i 4 and #24- This can be clearly seen 
in Fig. [21 where we plot the maximal and minimal values assumed by R e/I (left) and i? TM 
(right) as a function of E/E* in the absence of the muon energy loss due to synchrotron 
radiation, when varying the three-family mixing angles (#13, #23 and 612), the active-sterile 
mixing angles (#14, #24 and 6*34) and the CP- violating phases 81, 5 2 and 5 3 within the presently 



allowed region (see Ref. |35|). Notice that, when the parameters (a, Pi, (3%) are varied, these 
uncertainties give a little energy dependence to the observed ratio R efM but not to R T)1 , 
as it was expected from the previous discussion. The difference between the three- and 
four- family mixing is not large in R e/I , but it can be remarkable in R T)1 . 

In Fig. [3] we present the maximal and minimal values assumed by R efl (left) and R Tfl 
(right) as a function of E/E* in presence of the muon energy loss due to synchrotron radiation 



(taken from Fig. 15 of Ref. 22J)- In Fig.[3](left) we can see that, in presence of muon damping 
effects, the energy dependence of R efl becomes too large to disentangle new physics signals 
from the underlying uncertainties on the three-family mixing matrix parameters. Notice 



that theoretical uncertainties which are not taken into account in Ref. |22j can only worsen 
this result, making extremely difficult (if not impossible) to use R e/I to look for new physics. 
On the other hand, in Fig. EDright) we can see that the energy dependence induced in R TfJi 
by the muon energy loss is much milder than for R efl . This flavor ratio could still be used 
to distinguish three- from four-family mixing when #34 is near to its allowed upper bound. 
This is a further indication that R Tfl is the relevant observable to look for a signature of the 
(3+l)-scheme. 

Notice, however, that the energy domain for which muon damping effects are negligible 
is for E < E*. For E* ~ 10 5 GeV (taking e& = 1), the energy range for which it could be 
easier to use this flavor ratio to look for new physics is below the threshold for r production 
(E T > 10 6 GeV, see Refs. [lo|, Q). It should be stressed, on the other hand, that even for 
E > E* the three- and four-family R Tfl ratios do not overlap (for maximal #34) and therefore, 
with statistics large enough, the two models could be distinguished 3 . This is the main result 
of this paper. 

V. FLAVOR TAGGING AND STATISTICAL REQUIREMENTS 

It is now mandatory to discuss the statistical requirements needed to take advantage of 
high-energy v from astrophysical sources to look for signals of sterile neutrinos. 



This is not true for the R efi ratio. 
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FIG. 3: The same as Fig. [2] when the effect of the synchrotron losses of high energy muons are 
taken into account. The five cases are considered with the reference parameters = oo, 30, 3 in 
Fig. 15 of Ref. [^j for each case. 




First of all, we remind how different flavor could be identified in a neutrino telescope 
experiment. As it was shown in Ref. 14j, muon-like events can be identified with a good 
efficiency by looking for a single up-going muon entering the detector, whose energy can be 
estimated from the amount of Cerenkov light produced by photons and pairs radiated 
from the muon track. These are called " track events" . 

A second class of events are those in which the release of a large amount of Cerenkov 
light, generated by a shower induced by a neutrino interaction in the detector volume, is 
observed. This signal is produced mainly by v e {y e ) CC interactions, albeit with a (significant) 
contribution by and v T CC interactions and NC interactions of the three active neutrino 
flavors. A good knowledge of "track events" and "shower events" is needed to have a handle 
on the R eil flavor ratio (see, for example, Ref. 2l|). Moreover, the v e and z/ M fluxes must be 
known with a good precision to quantify the high-energy atmospheric neutrino background 4 
that is relevant for E v < 10 5 GeV. However, as it can be seen in Figs. [2](left) and [3](left), 
the ratio of "shower" to "track events" is strongly dependent on the spectral shape of the 
flavor fluxes. This introduces a further problem in the utilization of this observable to look 
for new physics beyond the three-family neutrino mixing scenario or to pin down the details 
of the astrophysical source of the high-energy neutrinos. 

Eventually, v T CC interactions can be singled out by looking for "double bang" or "lol- 
lipop" events [10j • The former are events in which two distinct energy releases are observed 



4 Notice that 7 x 10 5 atmospheric neutrino events with energies E v > 100 GeV are expected in 10 years of 
data taking at IceCube, [48 1. 
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inside the detector volume, with space and time separation consistent with the propagation 
and decay of an energetic r. The latter consist of events in which the second " double bang" 
shower is observed along the r track, but the first is missed. Both event classes are specific to 
v T CC interactions. Inverted lollipops (where the second shower is missed and the first one 
is observed) are not useful to identify r events, since they can be confused with an energetic 
Ufj, CC interaction in which a hadronic shower is reconstructed together with a muon track. 
The ratio of the sum of "double bang" and "lollipop events" to "track events" is, thus, a 
measure of R TlJL . 

We should now try to specify the statistical requirements needed to use flavor ratios at 
neutrino telescopes (measurable as explained above) to look for new physics from astrophys- 
ical high-energy neutrinos. 

Two sources have been considered in this paper: AGN's and GRB's. For both sources, we 
lack of an accurate model to predict the expected neutrino rates, and we can only compute 
their flavor ratios. This is useless, however, if we cannot establish if the observation of a 
deviation of flavor ratios from the three-family neutrino mixing expectation can indeed be 
a signal of new physics or it can just be ascribed to statistical fluctuations of the standard 
values. 

In Ref . [50(| , the present catalogues of TeV 7-ray sources were used to infer the expected 
neutrino fluxes to be measured at neutrino telescope experiments. Indeed, the hadronic 
interactions that are the source of astrophysical neutrino fluxes also create a large number 
of 7T° and i] particles. These particles decay eventually to photons, generating a high-energy 
7 flux. For a transparent source in which photons are produced mainly by this process, the 
relation that can be inferred on the neutrino flux by the observation of a high-energy photon 
flux is robust 5 . The neutrino flux is as large as the photon flux, at least. 

The three brightest 7-ray sources detected by the HESS telescope are the Crab neb- 



ula, RXJ1713. 7-3946 and Vela Junior [51], 152J. Their integrated flux above 1 TeV is 
(2.1,2.0,1.9) x 10~ n (cm 2 s)" 1 , respectively. The photon emission of the last two objects 
is believed to be of hadronic origin. Considering a photon flux half that of RXJ1713. 7-3946 
and Vela Junior, produced through hadronic mechanism, the corresponding neutrino flux 
above 1 TeV is 1 TeV) ~ 10~ n (cm 2 s) _1 . Using this neutrino flux, with an unbroken 

power law spectrum, the total number of e, \x and r contained events that can be observed in 
a km 3 water equivalent detector is 10.3, 9.6 and 2.9 per year, respectively. The expected flux 
of up-going i/-induced muons is <3> M = 5.6 (km 2 year) ~ 1 . These numbers take into account 



the different neutrino flavor energy spectra, as in Fig. 3 of Ref. [50j, the neutrino-nucleon 



5 It must be reminded, however, that photons can also be produced through inverse Compton scattering of 
relativistic electrons on radiation fields. In this case, the hadronic component is poorly measured and the 
neutrino flux can be much smaller than the photon flux. Eventually, for an opaque source, the photons 
can be absorbed inside the medium and subsequently emerge again at lower frequencies. 
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cross-sections and the Earth absorption, but not the detector efficiency to identify a given 



neutrino flavor. Similar results have been found in Refs. 53j, |54J |. 



The photon emission of the brightest AGN's present in the TeV 7-ray catalogue (Mkn421 



and Mkn501) 55j is strongly varying with time. In the emission peaks, their average photon 
flux above 1 TeV can be several times that of the Crab nebula. The extrapolation from 
the photon flux to the neutrino flux is difficult, however, since we have no indication on the 
hadronic or leptonic origin of the photon emission. 

Considering for AGN's and GRB's a flux similar to that of a galactic point source of 



the example given in Ref. [50], we should get 0(100) v e and events and O(30) v T events 



after 10 years of data taking at a km 3 detector. Similar statistics has been used in Ref. 56] 
to study neutrino decay at neutrino telescopes. It is easy to see that such statistics is 
too small to disentangle a (3+1) sterile neutrino effect from the statistical fluctuations of 
the three-family neutrino signals, even in the case of #34 near to its experimental limit, 
#34 ~ 35°. The uncertainties on the three-family parameters within the presently allowed 
region in R TI1 are 0(10%). The (gaussian) statistical fluctuations in R Tfl (after 10 years of 
data taking in the case of an ideal detector, given the numbers above) can be as large as 20% 
at la. The expected three-family R T/I flavor ratio is, thus, fdm ' ) ~ 0.30 ± 0.03(theo) ± 
0.06(stat), to be compared with i?^~ fam ' 1 (#34 = 35°) ~ c^r!^^^ ~ 0.2. This means that 
it is extremely difficult that (3+1) sterile neutrinos produced by a single source can be 
observed and unequivocally distinguished from theoretical and statistical fluctuations of the 
three-family mixing scenario in the next generation of neutrino telescopes. A factor 30 more 
statistics is needed to detect unambiguously a sterile neutrino signal with maximal 6*34. A 
possible way out of this problem is to integrate over all the galactic and extragalactic sources 
(or, for example, to sum over many GRB's events of similar intensities) 6 . In this case, it 
is possible that the next generation of neutrino telescopes could observe or constrain the 
four-family mixing model by looking at the R Tfl flavor ratio. No signal can be expected from 
the R efJi one. 

VI. CONCLUSIONS 

We have studied the effects of the (3+l)-scheme, which is constrained by all the experi- 
mental data except LSND, on the flux ratio of the high energy cosmic neutrinos. In principle, 
if a large number of events are accumulated in the future, there is still a chance to have a 
signature of the (3+l)-scheme. Furthermore, we have taken into account the theoretical un- 



certainties which were discussed in Ref. 22( in the case of high energy neutrinos from GRB. 
While these uncertainties make the prediction for the observed ratio R efl = (is e +i> e )/ (z^ + z^ t ) 
ambiguous, the ratio i? TM = {u T + z/ T )/(z/ At + z/ M ) is less sensitive to the uncertainties and is 



The authors thank P. Lipari for discussions which lead to these ideas. 
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sensitive to the sterile neutrino mixing angle 634. The energy spectrum of these flux ratios 
can also give us the information on the theoretical uncertainties, particularly that of the 
muon energy loss. While the statistical errors of data from a galactic point source in the 
next generation of neutrino telescopes are estimated to be too large to distinguish the three 
and four family schemes, if we can gain statistics by, e.g., summing over data from many 
sources, then it might be possible to have a signature for the (3+l)-scheme in observations 
of the high energy cosmic neutrinos. More detailed studies on the systematic errors to 
identify v T as well as the theoretical uncertainties would be necessary to derive quantitative 
conclusions. 



APPENDIX A: ANALYTIC EXPRESSIONS FOR THE FLAVOR FLUXES AND 
FLAVOR RATIOS AT SECOND ORDER 



At second order in 13 , 9u, 824 and 77 = 11/4 — 623, the picture is far more entangled. Both 
$14 and #24 are present, together with the three CP-violating phases. Notice that S13 always 
appears in combination with the phase factor cosy?, both in the expressions for the flavor 
fluxes and in the flavor ratios. We therefore define the reduced parameter S13 = S13 cos (p, 
where if = 5i — 5 2 + S3 as in eq. ffTTT) . However, at second order, also the phases 61 and 
5\ + 63 appear, and thus three independent CP-violating phases are present at this order in 
the flavor fluxes expressions. 

Notice that our results reduced to the corresponding three-family expressions at second 



order in 9 13 and 77 presented in Ref. |2l| for # 14 = # 24 = 6*34 = 



For the electron neutrino flux at the detector, we get: 
$ e /$° = {1 + (77 + s 13 cot 20 12 ) sin 2 20 12 



1 \[2 

2s 14 + 2 ( S 14 ~~ S 24) sin2 2^12 + ~Y Sil1 4#12-Sl4S24 COS \ ( A l) 



For the muon neutrino flux at the detector, we have: 



<V$° = jl - ~ (77 + s 13 cot 20 12 ) sin 2 29 12 

+ s 2 3 sin 2 26 l2 - r/sis sin 46 12 + rj 1 (4 - sin 2 20 12 ) 

1 \[2 
- 2s 24 - - (s 2 4 - s 24 ) sin 2 29i2 — sin 46>i 2 si 4 s 24 cos 81 )■ ( A2) 
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For the tau neutrino flux at the detector, we have: 



,. | ^ 1 - - (77 + S13 cot 20i2 ) sin 29 12 



sj 3 sin 2 26»i2 + VSis sm46 12 - 77 (4 - sin 2 20 i2 ) 
1 

S 24 _ 7 ( S 14 _ S Ia) Sin2 2 ^12 J" Sl n 4012SUS24 COS 5i 



77 (4 - sin 2 20i 2 ) - -S13 sin 40 12 



+ S 34 (2 S 14 + 3^34) 
+ (C34S34) |s24COS(5 3 

+ v2si 4 cos(5i + <5 3 ) sin 2 20 i2 (77 cot 20 i2 - si 3 )] 
Eventually, for the sterile neutrino flux at the detector, we get: 

$s = <M#34 -> #34 + 7T/2) 



(A3) 



(A4) 



Flavor ratios at first order in 6*13, #14, 2 4, 77 = vr/4 — 2 3 an d e(£' 1/ ) = 2 — X(E U ) have been 
presented in eqs. ( Tl5]) - (TTS1) . At second order in the same quantities we get: 



and 



<" fam) = 1 + \ (l - \ sin 2 20 i2 ) + ^(77 + 13 cot 20 12 ) 



sin 2 20 12 

2(0? 4 + 24 ) ( I - - - sin- 20 2 ) ^ sin 40 12 0i4024 cos 5x 



3 

X 



- ^0 2 3 sin 2 20 1 2(l + 3sin 2 20 1 2)+ 770i3 ^1 + |j sin 2 20 12 ) sin40 12 

- 4t7 2 fl-^sin 2 12 - ^sin 4 20 12 

3 

+ — e(7/ + 13 tan20 12 )sin 2 40 1 2 



^(4-fam) = c 2 



31 



I $//" ( 1 - i sin 2 20 12 J + i(er7 - 40 2 3 ) sin 2 20 12 



+ -(e + 8r/)0i 3 sin 40i 2 



+ sin 20 34 < ^p0i4 cos(5i + 5 3 ) (77 cos 20 12 — 0i 3 sin 20 12 ) sin 20 12 
v2 

277 M - ^sin 2 20 12 J - iesin 2 20 12 - ^0i 3 sin40 12 



+ 024 COS S 3 

+ {0 2 4 + 24(0 2 4 + 2 4 )} 



(A5) 



(A6) 
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